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Abstract 



We investigate the non-perturbative quantization of phantom and 
ghost degrees of freedom by relating their representations in definite 
^ ■ and indefinite inner product spaces. For a large class of potentials, 

we argue that the same physical information can be extracted from 
either representation. We provide a definition of the path integral for 
these theories, even in cases where the integrand may be exponen- 
tially unbounded, thereby removing some previous obstacles to their 
non-perturbative study. We apply our results to the study of ghost 
fields of Pauli-Villars and Lee- Wick type, and we show in the context 
of a toy model how to derive, from an exact non-perturbative path 
integral calculation, previously ad hoc prescriptions for Feynman dia- 
gram contour integrals in the presence of complex energies. We point 
out that the pole prescriptions obtained in ghost theories are opposite 
to what would have been expected if one had added conventional ie 
convergence factors in the path integral. 
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1 Introduction 



In this article we discuss the non-perturbative quantization of bosonic phan- 
tom degrees of freedom. A phantom is defined as a degree of freedom for 
which the sign of the kinetic term is opposite to that of an ordinary bosonic 
degree of freedom. Phantoms may often be quantized using either an ordinary 
Hilbert space representation or an indefinite inner product representation. In 
the latter case, the phantom is usually called a ghost. 

It is often physically required that we represent phantom degrees of free- 
dom on an indefinite inner product space. Examples include the timelike 
components of gauge fields, Faddeev-Popov ghosts, Pauli-Villars regulator 
ghosts, and the regulator ghost fields postulated in Lee- Wick [U ultravio- 
let completions of field theories. 

To describe such degrees of freedom non-perturbatively, Boulware and 
Gross [4j constructed a functional integral from an indefinite inner product 
representation of the canonical commutation relations. The resulting func- 
tional integral had two main problems, both of which we address in this 
article. 

The first problem occurs if the potential has odd terms. The integrand 
then becomes exponentially unbounded and the functional integral naively 
diverges, seemingly casting the non-perturbative existence of the theory in 
doubt p]. This problem already arises in the presence of source terms that 
are linear in the fields |5]. We show in this article that such linear - and 
in certain cases even quadratic - exponentially unbounded integrands can 
be treated non-perturbatively by defining path integrals in the distributional 
sense on appropriate test function spaces of a type studied by Gel'fand and 
Shilov in their theory of generalized functions [51 [7]. 

The second problem occurs when quantizing relativistic fields using the 
indefinite metric representation, since the region of integration may not be 
manifestly covariant. This happens for the electromagnetic field, as discussed 
by Arisue et al. |8j, and for the Dirac boson discussed in the present article. 
In particular, the region of integration for each ghost degree of freedom is 
effectively the imaginary axis. In general, there is no coordinate-independent 
way of distinguishing ordinary and ghost degrees of freedom, so that the func- 
tional integral is not manifestly covariant, though it may still be implicitly 
covariant [8]. The main result of the present article overcomes this problem 
by showing that information appropriate to the indefinite inner product rep- 
resentation can in fact be extracted in the usual way from the usual Hilbert 
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space representation, where covariance of the integrand is manifest. 

For gauge fields, the existence of a well-behaved analytic continuation to 
Euclidean time makes these issues less urgent. However, for fields such as the 
Dirac boson, which are of particular interest as Pauli-Villars and Lee- Wick 
ghosts, the Euclidean functional integrand is exponentially unbounded. For 
these fields, the real-time formalism works fine, as we illustrate. Astonish- 
ingly, for fields that appear at most quadratically in the action, even the 
exponentially unbounded Euclidean functional integral can be calculated in 
an appropriate distributional sense, as we discuss in section [51 

It is in fact common practice to derive results based on heuristic functional 
integrals, in theories containing indefinite-metric gauge field components or 
ghosts, by using the classical covariant action. One may justifiably ask why 
not. In the real-time formalism, the integrand is exp(i5'). As long as S is 
real, negative terms in 5* are not necessarily worse than positive terms in S, 
and the oscillatory integral may be equally well defined in either case. 

In this article we justify the use of the ordinary covariant action and 
integration region for potentials sufficiently general to cover the applications 
mentioned in the introductory paragraphs. In particular, we show that the 
indefinite-metric ground state expectation values can be obtained from the 
definite-metric transition function using the usual T —ioo prescription. 

Because they lacked a non-perturbative definition of the indefinite inner 
product theories corresponding to the unbounded path integrands that we 
are now able to treat, various authors [H El E] studied the consistency of ad 
hoc prescriptions for defining the propagators, and more complex diagrams, 
order by order in perturbation theory. Lee and Wick [H [2] introduced one 
such prescription, requiring the contour integrations defining various dia- 
grams to be continuously deformed to avoid the movement of poles in the 
complex plane as one changes the parameters of the theory. Cutkosky, Land- 
shoff, Olive and Polkinghorne ^ attempted to generalize this prescription to 
coalescing singularities not covered by Lee and Wick. 

The framework of this article allows a non-perturbative functional integral 
to be calculated in the distributional sense for various kinds of unbounded 
path integrands, including the ones studied by Lee and Wick. In principle, we 
expect the ad hoc prescriptions of the above authors to be either reproduced 
or corrected in our non-perturbative framework. We motivate this statement 
with a couple of simple exact calculations in a toy model, where we succeed in 
deriving Lee- Wick-type pole prescriptions from first principles via explicitly 
non-perturbative calculations. 



3 



2 Related work and applications 



In addition to the seminal references discussed in the introduction, we also 
point out the following related work. 

The study of phantoms has recently been of interest in cosmological model 
building. They are used, for example, in some models of cosmological dark 
energy. The literature on the subject is too large to cite representatively, but 
see for example [9J. There are open questions regarding the consistency and 
stability of these models [101 E] ■ It is our hope that the non-perturbative 
methods introduced in the current article may throw further light on the 
subject. For example, in section [S] we learn how to non-perturbatively treat 
unbounded and complex potentials, which may be useful in studying, from 
a quantum-mechanical point of view, the stability questions discussed in, for 
example, [12] and [T3] . 

In [T^ [15] Erdem, and in [16] 't Hooft and Nobbenhuis discuss a novel 
kind of symmetry transformation consisting of a rotation of real positional 
coordinates to the imaginary axis, with the aim of ruling out a cosmolog- 
ical constant. The rotated representation that these authors use for their 
non-relativistic particle toy modle is identical to the one we discuss in sec- 
tion [3] for the indefinite inner product theory. These authors are particularly 
interested in the relationship between the real and imaginary coordinate rep- 
resentations. Since we study this relationship in detail in the present article, 
it is conceivable that our mathematical framework may have further appli- 
cations in this direction. 

Also related to the cosmological constant problem is the paper [17], which 
introduces phantom fields to cancel the ordinary matter contribution to the 
vacuum energy. Again, our non-perturbative approach to phantom fields 
may be useful in the study of these models. 

In the paper [18], Antoniadis et al. promote the study of theories with 
ghosts in the real-time formalism, as opposed to the Euclidean formalism. In 
the real-time formalism, they study the issue of ie prescriptions for resolving 
ambiguities in Feynman diagrams. However, we disagree with the assump- 
tion of their approach, which derives pole prescriptions by introducing ie 
terms to provide convergence factors in path integrals, and states (we believe 
incorrectly) that the convergence terms are necessary and that theories with- 
out such convergence terms are ambiguous. As we discuss in section [9l the 
pole prescription for ground state expectation values in ghost theories is in 
fact opposite to what one would expect from such a convergence factor. Our 
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derivation is unambiguous, and does not require a convergence factor, since 
our path integral is well-defined without it. As a result of their assumptions, 
the authors of [T^ find opposite iO prescriptions for the two-point functions 
of a particle and a corresponding ghost, in confiict with our results, which 
agree with many authors starting with Pauli and Villars |19J. In the light of 
these remarks, we believe that the conclusions of [181 ^lay have to be revis- 
ited or reinterpreted. In this regard, we believe that the expectation values 
calculated in [18] may be valid in a representation based on a vacuum that 
is not a ground state. It is unclear to us how such a vacuum can be defined 
invariantly. 

In section [T21 we provide further non-perturbative calculations demon- 
strating what happens when the poles may travel in the complex plane. 

There is a large body of work in scattering theory that includes res- 
onances in the description of quantum systems via complex-coordinate or 
complex-momentum methods that displace either the position or the mo- 
mentum representation into the complex plane [20l [21], [221 l23l [211 [25]. In 
these cases, the original and displaced representations are not equivalent, yet 
they are expected to encode the same physical information. What we do 
in the present article may be seen as a special case of this method, made 
rigorous for a specific class of potentials. We show that the same informa- 
tion is encoded by the representation based on the definite inner product, 
whose configuration space is the real line, and the representation based on 
an indefinite inner product whose configuration space is the imaginary axis. 

In scattering theory, so-called Siegert or Gamow states may be used to 
represent resonances. These are states with complex momentum [20], which 
may be given precise mathematical meaning in the framework of section 
\5\ of the current article, where such states are defined as distributions on 
test function spaces of Gel'fand-Shilov type. Since we are mainly interested 
in field theory applications where interactions are polynomial, we do not 
treat sufficiently general potentials for our results to be directly applicable 
to many traditional non-relativistic scattering problems, but perhaps the 
mathematical machinery can be generalized. 

The author of [2H] describes a treatment of these states in the context of 
rigged Hilbert spaces. His construction should be very closely related to the 
one of the current article. 

We also mention the following recent work on alternative physical inter- 
pretations of theories of Lee- Wick type, which differs from the approach of 
the original authors [27]. 
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3 Relating definite and indefinite quantiza- 
tions: The free particle 

We will show that certain quantum systems can be quantized in either a 
definite or an indefinite inner product representation. We will also argue 
that, under certain conditions, the choice is a matter of convenience, since 
physical quantities appropriate to one representation can be obtained from 
results calculated in the other. 

To introduce the ideas, in this section we start with the ordinary free 
particle in one dimension with Hamiltonian 

n2 



2m 



The transition function is given by 



^ik{y-x)-itk^ /2m. 

2tt 



V 27Tt 

where the last line is calculated for positive t. Here we have summed over a 
complete generalized basis 

(x|fc) = e*'=^ p\k) = k\k), {k'\k) = 27T5{k' -k), keK, 

of a Hilbert space with positive definite inner product. The energy spectrum 
is positive, and the oscillatory integral has a rigorous interpretation as the 
Fourier transform of a distribution. 

Notice that the form ([T]) of the spectral representation may be continued 
to a function analytic in t on the lower half of the complex plane. 

The same Hamiltonian can be quantized in an indefinite inner product 
representation. We will recast the representation used by Boulware and Gross 
[5] in a more convenient notation similar to the one used by Arisue et al [S] . 
This proceeds by representing q and p as Hermitian operators in an indefinite 
inner product space spanned by generalized bases 

p\ik)) =ik\ik)) , {{-ik\ik')) = 2n 6{k' - k), keR, 
q\ix)) = ix\ix)) , {{—ix\ix'))=6{x' — x), x G R, 
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where 



{{ix\ik)) = e-^*^^ 

Note that the eigenvalues of the Hermitiar0 operators q and p occur in com- 
plex conjugate pairs, which is allowed in an indefinite inner product space 
[28] . The corresponding pairs of generalized eigenstates are dual null states. 
The configuration space is now the imaginary axis, where wave functions take 
the values 

(f)(ix) = {{ix\(t>)) . 

Thus, a generic wave function in the position representation has as its do- 
main the imaginary axis. There is no assumption of analyticity - we cannot 
analytically extend a generic wave function away from the imaginary axis, 
so that {{z\(j))) is in general undefined unless z e iR. We will, however, have 
reason below to investigate subspaces consisting of test functions that are 
analytic. 

The completeness relations in terms of these sets of states are 

f dk 

dx I —ix)) {{ix \ = 1 = / — I —ik)) {{ik \ 



27r 

corresponding to the inner product 

dx {ijj{—ix))* (j){ix). 
satisfying 

(W)) = ((«)*• 

This inner product is indefinite. The state space is not a Hilbert space 
but a Krem space, which may be constructed by completing the space of 
square-integrable functions on the imaginary axis with respect to an auxiliary 
inner product to obtain a topologically complete vector space, which 
is then equipped with the above indefinite inner product, discarding the 
auxiliary CP' inner product [2S]. We emphasize that all that matters about 
the auxiliary £^ inner product used in the intermediate step is the resulting 
topological vector space. The auxiliary inner product is not unique and has 



^Hermitian operators on indefinite inner product spaces are sometimes also called 
pseudo-Hermitian. We will not use this qualifier. 
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no physical meaning. The construction is analogous to using an arbitrary 
and non-physical Euclidean auxiliary inner product to define the topology of 
Minkowski space. 

Our goal is to investigate the relationship between definite and indefi- 
nite quantizations of the same Hamiltonian. We will therefore compare the 
ordinary Hilbert space transition function {y\ e~*^* \x = 0) to the quantity 
{{y\ e~*^* |x = 0)) calculated in the indefinite representation, where a suit- 
able meaning for ((?/|, y G R, will be assigned below. 

At first glance, this might not seem like a reasonable thing to do. The 
two quantities are distributions on different test spaces. As distributions, 
they are therefore incomparable. However, as long as there exists a range of 
t for which the distributions are representable by continuous functions in y, 
we can compare these functions. 

The first observation is that we need a final state {{y\, y G R, whereas 
generic states in the indefinite representation are functions (f){ix) = {{ix\(j))) 
defined on the imaginary axis. Generalized functions such as {{y\ can, how- 
ever, be defined provided we can find a subspace of the state space consisting 
of test functions that may be analytically continued from the imaginary axis 
into the complex plane and are, at the same time, invariant under time evo- 
lution. For the free particle, a simple test function space satisfying these 
conditions is the Gel'fand-Shilov space Sy'^i'^^) [6], to be discussed in more 
detail in section \5[ This test space consists of entire functions satisfying the 
growth conditioio 

\y''(P{x + ty)\< Ck exp {-a\y\^ + b\x\^) , (2) 
where Ck, a, and b are positive real numbers that are allowed to depend on 

■^The space S^"^ defined by Gel'fand and Shilov is based on imposing growth conditions 
on the real axis and corresponds, in the general notation introduced in the section [51 to 
Sljlix). The space S^j'^{ix) used here is simply obtained by imposing the same growth 
conditions on the imaginary axis, exchanging the roles of x and y. 

•^Although the powers of x and y in the growth conditions are preserved under time 

1 /2 

evolution, it should be noted that the sign of b is not, so that the space S^^^ {ix) is not quite 
closed under e~*^*. This makes it necessary to consider instead families of generalized 
functions on subspaces denoted by 5*^^2^(13;) that are mapped into each other under time 
evolution, where A and B are related to a and 6 in (U) . The method is developed rigorously 
by Gel'fand and Shilov in [7j, and can be applied as long as the powers of x and y in the 
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For arbitrary z G C, we may now define the generahzed function {{z \ as 

((^10)) <peS\'/l{tx). 

1 /2 

Acting on an element of Sy^{ix)^ we have, purely formally, 

((y|e-^*|0))= ! dk{{y\e-^'''\-ik)){{im) 



dk{{y\e-''''\-ik)) J dx{{ik\-ix)) {{ix 

However, neither of the supposed kets | —ik)) nor | —ix)) is a test function, 
and the above should be read as a shorthand for the rigorously defined 



nentially growing factor e^^ , the integral over k converges, since S\^l^{ix) is 
closed under the Fourier transform [6], schematically J^{S\^l^{ix)) = Sljl{ik), 



((l/|e-'^*|0)) = / dke"' '+''y{{tk\^)) 

It is here that the utility of Sl^^l{ix) comes to the fore. Despite the expo- 

nfofTT-ci r^^rQY k couvcrgcs, since S-^^^ 
latically JF(S'^^2(^2;)) = 
which means that {{ik\(j))) decreases as e"'^^^ for some positive h. Then 

where the added convergence factor, entirely superfluous in the convergent 
integral on the first line, has allowed us to exchange integrations in the second 
line. 

We see from the second line that the generalized function ( {y \ e~*^* can 



growth condition are preserved under time evolution. In the interest of brevity, we will 
commit a slight abuse by leaving this complication implicit in the rest of this paper. 
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be represented by the kernel 



{{y\e'''''\-ix))=\im / dk {{y\e-'''' \-ik)) e-'" {{ik\-ix)) 

€— »0 



lim j dke'^^'+^^y-'''^-'^\ (4) 



= e^'^/^W — e*'"^^""'^''/^*, (5) 
V 27rt ' 

where again the formally undefined inner products on the first line should 
be taken as a useful mnemonic for the well-defined expression on the second 
line, copied from ([3]). Again, the last line was calculated for positive t. 

We have obtained a spectral representation of the kernel {{y \ e~*^* | —ix)) 
in the indefinite inner product representation. The expression (jlj) in fact 
defines a function analytic in t on the upper half plane. 

For X = 0, we may now compare this with the definite representation, 
and we find 

((y|e-^^*|x = 0)) =z (yl e-'^* |x = 0) 
for positive t. For negative t, one may check that 

((y|e-^^*|x = 0)) = -z (y|e-^* |x = 0). 
To summarize, we have obtained the following structure: 

• The function i {y\ e"*-'^* \x = 0), calculated in the definite inner prod- 
uct representation, has a positive-energy spectral representation that 
defines a function analytic in t on the lower half plane. 

• The function {{y \ e~^^^ \ x = 0)), calculated in the indefinite inner prod- 
uct representation, has a negative-energy spectral representation that 
defines a function analytic in t on the upper half plane. 

• On the positive real axis in t, these two functions coincide. Together, 
they therefore define a single function analytic on the entire complex 
plane except for a cut on the negative real axis, where the value of 
{{y\ e~*^* |x = 0)) is obtained by approaching the cut from above, and 
i {y \ e"^^^ |x = 0) is obtained by approaching the cut from below. 
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We see that, for positive real t, the results calculated using the two representa- 
tions coincide]^ If we knew beforehand that the two spectral representations 
determined the two halves of a single function analytic in t, we could have 
used either representation to infer the result for the other. In section [5], we 
will discuss under what conditions this can be done. 



4 The phantom free particle 

We briefly discuss the quantization of the particle with opposite sign Hamil- 
tonian 

2m 

Now the positive-definite inner product momentum eigenstates have negative 
energy 

.2 



and 



e 



^ ^ik{y-x)+itk^ /2m ^g^ 

27r 



V 27lt 

Note that, although the energy is unbounded below, the real-time represen- 
tation exists. The spectral representation ([6]) now directly defines an analytic 
function in t on the upper half plane. 

The same transition function can be written in terms of a positive energy 
spectral representation using the indefinite inner product space defined in 
section [31 The states 

\ik)) , keK, 



Since the inner product in the indefinite representation has no canonical choice of 
overall sign, we could have made these two quantities coincide for negative i, and positioned 
the cut at positive t, by defining the inner product with an opposite sign. 
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are now eigenstates of H with positive energy k'^/2m. Similarly to section [31 
we find the spectral representation 



and ([7]) extends analytically to the lower half plane. 
We have found the structure: 

• The function —i {y\ e""^^^ \x = 0), calculated in the definite inner prod- 
uct representation, has a negative-energy spectral representation that 
defines a function analytic in t on the upper half plane. 

• The function {{y \ e"*-'^* \ x = 0)), calculated in the indefinite inner prod- 
uct representation, has a positive-energy spectral representation that 
defines a function analytic in t on the upper half plane. 

• On the positive real axis in t, these two functions coincide. Together, 
they therefore define a single function analytic on the entire complex 
plane except for a cut on the negative real axis, where the value of 
{{y\ e~*^* |a; = 0)) is obtained by approaching the cut from below, and 
—i {y\ e~*^* \x = 0) is obtained by approaching the cut from above. 

Again, for positive real t, the results calculated using the two representations 



5 General potentials 

We now generalize the above analysis to a wider class of potentials. We shall 
see that the transition functions in the definite and indefinite quantizations 
typically describe a single analytic function for different regions of complex 
time that overlap for an initial interval of real time. In sectional we will apply 
this correspondence to the extraction of ground state correlation functions. 

The general problem will be as follows. Consider a one-dimensional 
Hamiltonian 




ix)) 




coincide. 
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that is Hermitian with respect to the ordinary Hilbert space inner product 
on the real hne, which requires that /i be real and that V{x) be a real- valued 
function of x. Our discussion will include the case that is of interest in 
phantom and ghost theories, where the spectrum of H, which is real in the 
ordinary Hilbert space representation, may be unbounded below. As a result, 
the usual rotation to Euclidean time may not exist |f| but typically there is 
no obstacle to calculating the real-time transition function, which will be 
denoted by 

(y|e-^*|a;). (8) 

In the indefinite inner product representation constructed in section [HI where 
the configuration space of q is the imaginary axis, the operators p and q 
are still Hermitian, as is the Hamiltonian. However, the spectrum of the 
Hamiltonian may be complex]^ Of most interest will be those phantom or 
ghost theories for which, in the indefinite representation, the real part of the 
spectrum becomes bounded belowJil We will assume that the imaginary part 
of the spectrum, if any, is sufficiently well-behaved that e~*^* exists for all 
real t. The transition function in the indefinite representation is denoted by 

{{^y\e-^'''\-^x)). (9) 

Due to the different configuration spaces, the quantities ([8]) and ([9]) are not di- 
rectly comparable. For a comparison, we need states of the form {{y \ instead 
of {{iy\, whereas general states in the indefinite representation are functions 
(l){iy) = {{iy\(f))) that are only defined on the imaginary axis. Generalized 
functions such as {{y\ can, however, be defined on suitable test function 
spaces that are invariant under time evolution and consist of functions that 
may be analytically continued away from the imaginary axis into the complex 
plane. The choice of test function space is dependent on the Hamiltonian. 

Gel'fand and Shilov [HIE] introduced families S^{x) of test function spaces 
on the real line that are suitable for our purposes. These are defined, for 

^Although wc shall see that the relevant unbounded Euclidean path integral may in 
fact be calculable as a Gel'fand- Shilov generalized function. 

^We remind the reader Hermitian operators in indefinite inner product spaces may have 
complex eigenvalues. These always occur in complex conjugate pairs. The corresponding 
eigenstates are dual null states [2H]- 

^More generally, our arguments remain valid even in cases where the real part of the 
spectrum is unbounded below, as long as the density of generalized eigenvalues decreases 
sufficiently rapidly in the negative real direction. 
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01-, 13 < 0, as consisting of all infinitely differentiable functions on the 
real line satisfying the growth conditions 

where A, B and C may depend on cj). Writing the independent variable x 
in parentheses is not part of the original notation but will be necessary for 
later disambiguation when we discuss similar spaces based on functions 
defined on the imaginary axis. There exists an interpretation of the limit of 
infinite a or /?, and it is conventional to omit either index when it is infinite. 
The familiar Schwartz space is the same as = S'^{x). For other values of 
a and /?, the space is more restricted than Schwartz space since the growth 
conditions are stricter, and in fact for j3 < 1 these conditions are so strict 
that the functions (j){x) can be extended to analytic functions (f){x + iy) on 
some complex domain. 

The domain depends on (3. Specifically, for (3 = 1, the test functions (j){x) 
may be analytically extended to analytic functions on a strip 

{x + iy\x G R, \y\ < B}, 

where the width B of the strip may depend on 0. For f3 < 1, the test 
functions can in fact be continued an entire analytic function on the whole 
complex plane. In this case, the spaces S^{x) are characterized completely 
by the growth conditions 

\x''(f){x + iy)\< Ck exp (-a|a;|^/" + febr/^^"^^) , (10) 

where Ck, a, and b are positive real numbers that are allowed to depend 
on 0. Smaller a or (3 correspond to smaller spaces of test functions, and 
therefore larger dual spaces of generalized functions. The restrictions become 
so stringent that the test spaces are trivial for 

a + (3 <1. 

Under the Fourier transform, it may be shown that [6] 

J'{S^ix))=S-,ik). (11) 

In the indefinite representation, the wave functions are not defined on the 
real line but rather on the imaginary axis, and we will denote the spaces 
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where the same growth conditions are imposed in the imaginary direction 
by S^iix). As above, for /3 < 1 these test functions can be continued away 
from the imaginary axis to obtain entire functions. For (3 = 1, the strip of 
analyticity will now be a vertical one that includes the whole imaginary axis. 

The relation (|TT1) remains valid if we define the Fourier transform on the 
imaginary axis as 

T[(t)\{ik) = I dxe-'^''<P{ix). 



We will restrict our attention to potentials for which a test space can be 
selected on which the Trotter product formula 

e-*^* = lim { e--fV2-+^.MPe--v\^ , ^ 1 

N^oo\ J N 

holds at least for a finite time interval. This formula can be taken a basis for 
the construction of a path integral representation. If it holds, a generalized 
path integral may be defined by writing the (generalized) transition function 

((yl e~^^^ as 



If one writes the Fourier transforms, each acting to the left on a generalized 
function, as formal integrals, this can be seen to coincide with the introduc- 
tory textbook definition of the path integraljfl For related definitions of path 
integrals in terms of transforms, see for example [29] and the review in |30j . 

This expression will exist as a generalized function, but may not have a 
representation as an ordinary kernel. In other words 

((y|e-^*|-zx)) 

may not exist as an ordinary function of y and ix. It is only when the kernel 
can be represented by a function continuous in y and in a; at x = that we 
can compare 

((y|e-^*|-2x)) 



®It should be remarked that even the simplest textbook real-time path integral calcu- 
lations give meaning to the oscillatory integrals involved in a way that is essentially equiv- 
alent to calculating Fourier transforms of distributions as we do here. By using smaller 
invariant test spaces, our spaces of distributions are larger, allowing us to calculate a larger 
class of path integrals. 



15 



and 

(y|e-^* \x), 

in the point x = 0. Otherwise, as distributions, the objects 
and 

are incomparable, since they act on different test spaces. 

For the above construction, we saw that it is necessary to consider test 
spaces consisting of analytic functions that are invariant under all of e~*'^^^'^\ 

The appropriate test spaces will depend on the potential. We will choose 
suitable candidates from the family of spaces S^{ix), whose elements are 
either entire functions for /3 < 1, or analytic functions in a vertical strip 
including the imaginary axis for f3 = 1. 

We start with the spaces of entire functions S^{ix) for /? < 1. These 
spaces will be suitable only for potentials such that e"*''^^*^) can be extended 
to an entire function of z, since the space must be closed under this multiplier. 

The indices a and (3 are then further restricted by the growth properties 
of e~"^(^\ We will consider the class of potentials satisfying the inequality 

|e-^(^)| <Ce''l^^ (13) 

on the complex plane, which is true in particular if \^ is a polynomial of order 
p. Then a sufficient condition for S^ix) to be closed under this multiplier is 
that p < 1/a andp < 1/(1 -(3), oB 

«<-, /?>^. (14) 
p p 

Next, consider invariance of these spaces S^{ix), still for (3 < 1, under 
gitiep /2m+i/ip^ which is equivalent to invariance of the Fourier transformed 

^When the inequality is saturated, the space is not necessarily closed under the 
action of the multiplier. In this case, the appropriate test function spaces are families of 
subspacesS*^'^ of S^. As noted in a previous footnote, this technical complication does 
not affect the general argument and will be glossed over in this article. For full details on 
the technique, see [7]. 



16 



space Sp{ik) under e±*f(*'=)^/2m+ieM(«fc)_ \Ye will restrict attention to the fam- 
ily (I14p selected by the growth condition (IT^ . First, if the order p of the 
potential satisfies p > 1, then the first inequality constrains a < 1/p < 1, 
so that the Fourier-transformed space also consists of entire functions. This 
space will be invariant if 

(^<\. «>^- (15) 

Solutions to both sets of conditions f|T^ and f|T5|) may be found for potentials 
of order 1 < p <2. We may choose any such a and (3 satisfying the additional 
condition «+/?>! necessary for non-triviality of the test space. The smallest 
test space satisfying all these conditions for all p in this range is given by 
S\^l'^{ix)^ which happens to also have the nice property of being invariant 

1 /2 

under the Fourier transform. In fact, it is easily seen that S^i^{ix) satisfies 
the conditions for the extended range < p < 2, so that we may use as a 
test space for the whole range 

S\'jl{ix), 0<p<2. 

In fact, the invariance conditions are unchanged for arbitrary complex coef- 
ficients of p, p"^, X and x"^ in the Hamiltonian, enabling us to calculate path 

1 /2 

integrals as Gel'fand-Shilov distributions on S^^^ (^^) even when the integrand 
may be exponentially unbounded. 

We note that fll4p and (IT^ cannot be satisfied for polynomial potentials of 
order p > 2. As a result, none of the Gel'fand-Shilov spaces of entire functions 
may be used for these potentials. However, for our later arguments, it will 
be good enough to have test functions that are analytic in a vertical strip 
containing the imaginary axis, a property satisfied by the spaces S^{ix). This 
family of spaces, with suitable restrictions on a, will prove to be suitable for 
even potentials of order p > 2. 

According to [6], a sufficient condition for a function f{ix), defined on 
the imaginary axis, to be a multiplier with respect to a specific family of 
subspaces S^'^{ix) of S^{ix), is that it satisfy the estimate 

<Ce5'^g^''e^'"''^^ (16) 

for any e > 0. Consider, then, the term V{q) = cq^ in a polynomial potential, 
where c is real. Its contribution to the multiplier e^^^"^ is 



17 



We see that p must be even to satisfy the sufficient condition above. By 
writing down the ffist few derivatives, it is easily seen that 

\f\ix) < + |xn. 

This is not the tightest possible estimate, but it will suffice for our purposes. 
Using the Stirling estimate for large q 

In (g!) ~ g In g — g, 

we find 

\P{ix)\ < C q'^ e-'' p'' {1 + \xD 

Comparing with the estimate (fT6l) . we read off that any 

/3>1, a>0, 

will do, for any even power p in the potential. Of these, only (3 = 1 gives 
a space of analytic functions. Applying the same argument to the multi- 
plier e^*^*^*^-* on the Fourier transformed space gives the (sufficient) 
condition 

p>0, a>l. 

The smallest space of analytic functions satisfying both these conditions is 

which, again, has the nice property of invariance under the Fourier transform. 

In fact, while the above conditions are sufficient, it may easily be seen 
that the linear potential or source term e*'*'^ is also a multiplier on families of 
subspaces of Sl{ix), despite its exponential growth in the imaginary direction 
if A is real, so that we may allow also source terms in our potential. This is 
not true for higher odd powers. 

In fact, it may be checked that both test spaces that we have chosen are 
invariant under general source terms of the form 

fip + Ag, 

where /i and A may be complex. 

To summarize, we find that the following test spaces are suitable: 
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• H = up^ /2m + jip + V{q), where /i, z/ and V{q) may be complex, as long 
as |e~*''^*^^)| < Ce^l^l*' for < p < 2: A suitable test space is S'^y|(ix). 
Elements of the test space are entire functions. 

• H = /2'm+ fip+ Xq+V (q) , where fi and A may be complex and V{q) 
is an arbitrary even real polynomial: A suitable test space is Sl{ix). 
Elements of the test space are analytic on a vertical strip including the 
imaginary axis. 

A useful feature of the current formalism is that it allows the study of source 
terms in phantom or ghost theories. These caused some difficulty in prior 
treatments [HE]. In particular, while the exponential of a source term e*'^^ 
blows up on the imaginary axis where the states of the ghost theory are 

1 /2 

defined, it is still a legitimate multiplier on S-^^^^^ix), where its exponential 
growth is dominated by the faster rate of decrease of the test functions, and 
on families of subspaces of Sl{ix), where it just modifies the coefficient, and 
not the exponential power, of decrease of test functions. 

In both cases, the path integral may be defined, at least for a finite initial 
time interval, via the formula / fllj) . which defines a legitimate distribution 
despite the fact that the integrand in the naive path integral may be exponen- 
tially unbounded. 

A similar analysis can be made for the definite inner product version 
whose configuration space is the real line by replacing ix by x in the argu- 
ments. 

Potentials for which both Slj'^{ix) and Sl{ix) are suitable include the 
free particle and the harmonic oscillator, as well as the linear potential, of 
interest in some cosmological models, and the inverted harmonic potential 
that appears in two-dimensional string theory, despite the fact that the po- 
tentials in the latter two cases are unbounded below. The path integral may 
be defined for these potentials via the formula f[T^ . 

While Sl{ix) allows the treatment of higher polynomial potentials that 
cannot be treated using Sy^i^^)} we note that Sy'^{ix) allows us to treat 
some potentials that cannot be treated using Sl{ix). In particular, while 
Sl{ix) is good for real polynomials, the growth condition le"**^^*^^-*! < C e^'^'^ 

1 /2 

for S-^^^^{ix) does not require the momentum terms or V{q) to be real. 

In this case, the expression [W\] defines a good path integral even though 
the integrand may exponentially unbounded with growth of order e"'^ . 
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Given a test space consisting of analytic functions, a generalized function 
{{z \ can be defined as 

((#)) = 

for any z within the common domain of analyticity of the test functions 0, 
even though the full space of states is defined only on the imaginary axis. 

We will consider the phantom Hamiltonian introduced at the beginning 
of the section 

2 

H^-^ + ^^p + V{q). (17) 

and argue that under certain assumptions the quantities 

(y| e"*^* \x = 0) 

and 

^((l/|e-^^*|x = 0)), 

coincide for some range of y and for some initial time interval < t < T. 
This will hold for any real y when we use the test space Sy^i"^^)^ whose 
elements are entire functions, or for ?/ in a neighbourhood of in the case of 
Sl{ix), whose elements are analytic only on a strip including the imaginary 
axis. 

Our first assumption is that the distribution {{z \ e"*-'^* can be represented 
by a kernel, denoted by 

{{z\e-''''\-ix)) , 

that is a differentiable function in ^{z), '^{z) and x and a differentiable func- 
tion of t for some initial time interval Given this assumption, 
we may write 

((.|e-'"|^)>^/..((.|e-'«|--))^(-.0) 

for test functions (f) G Sljl{ix) or G Sl{ix). Here (t){z,t) is analytic, since 
the space of test functions was chosen to be closed under time evolution. 
Differentiating on both sides with respect to dz, and using the fact that the 
spaces S^iix) are "sufficiently rich" [6], we conclude that 

d-z{{z\e-'^'\-ix)) = 0. 

-"^^The example of the harmonic oscillator, for which the transition function rcfocuses to 
a delta distribution after a half period, shows the need for the time interval restriction. 
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In other words, {{z\e | ix)) is analytic in z. 

Due to this analyticity of {{z \ e"*-'^* the differential equation 

{idt-H {-dy, ly)) {{ly I e-^^* | x = 0)) = 0, 

which is satisfied by assumption for < t < T by the transition function in 
the indefinite representation (g ix, p ^ —dx) on the imaginary axis, is in 
fact true on the whole domain of analyticity in z. We find 

_ H {-id,, z)) {{z I e-*^* I X = 0)) = 0, 

1 /2 

on the entire real axis in the case of Sy^ iw) ! some real neighbourhood of 
zero in the case of Sl{iy). But for real values of z, this just coincides with the 
ordinary Schrodinger equation giving the time evolution of {y \ e~^^^ |x = 0) 
in the ordinary Hilbert space representation. 

Therefore, on the interval < t < T, the objects 

z((y|e-''^*|a; = 0)) 

and 

are ordinary functions satisfying the the same differential equation in y. They 
would therefore coincide if they satisfied the same initial conditions. However, 
the initial conditions at t = in the two representations are not ordinary 
functions, but rather delta distributions on two different test spaces. These 
distributions are not comparable. 

We therefore recast the problem in a form that gives an ordinary function 
as initial condition. This is done by considering the quotient 

^ {{z\e-^"'\x = Q)) 
^ ' ' ~ ((z|e-^ot|a; = o))' 

where Hq is the free phantom Hamiltonian considered in section HI and the 
denominator is the analytic continuation to all 2; G C 



{{z\e~'''^'\x = Q)) = e-'^'\ —e-'^^'^'^K (18) 



of the transition function 



((zi/|e-^^«*|x = 0)) = e-^-/%/£^e^'"^'/2*^5(l/) as t 0, 
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of the phantom particle on its imaginary axis configuration space. Given the 
form fll7p of the Hamiltonian, the differential equation satisfied by G{z, t) is 



{itdt — tH {—idz, z) + izdz — fxruz) G{z, t) = 0, 



(19) 



for < t < T, with initial condition the ordinary function 



G{iy,t)^l, t^O, 



on the imaginary axis, assuming that the effect of the interaction becomes 
negligible in the limit as t — ^ 0, so that the behaviour of the numerator 
approaches that of the denominator in this limit. For the class of smooth 
potentials under discussion, we expect this to hold, although we do not have 
general proofl"] 

To compare this solution with the one obtained in the definite representa- 
tion, we need to know what happens to G{z, t) for real values of ^ as t — > 0. 
To this end, we note that the solution of this differential equation gives a 
family of functions, indexed by t, that are analytic on an open domain in the 
complex plane in z including either a vertical strip containing the imaginary 
axis or the whole complex plane, depending on the choice of test space. This 
family approaches to the constant function 1 on the imaginary axis. If this 
family is uniformly bounded on some open neighbourhood of the origin for 
small enough t, we may use Vitali's theorem |3T[ |32] to conclude that the 
family converges to the constant function 1 not only on the imaginary axis 
but on this whole neighbourhood and, by implication, also on the segment 
of the real axis included in the neighbourhood. 

For the class of potentials under discussion, we expect the condition of 
uniform boundedness to be true, although we have not found a proof. We 
will in the following assume this assertion without proof. 

We would like to compare G{y,t) to the corresponding quantity 



in the definite representation, which satisfies the differential equation 



It is physically reasonable that the potential requires time to make its presence felt, 
so that the behaviour approaches that of a free particle in the zero time limit. 



G{y,t) 



{x\ e \x = 0) 
(?/|e-^^o* |a; = 0)' 



{itdt - tH {-idy, y) + iydy - jimy) G{y, t) = 0, 



(20) 
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for < t < T, with the initial condition 

Since the equations f|T9l) and fl20l) coincide for real z and have the same 
initial condition on the real line, their solutions coincide on the indicated 
time range. For this range, we conclude 

G{y,t) = G{y,t), (21) 

From their explicit expressions, the free transition functions in the two rep- 
resentations can be compared to obtain 

(d/l e-'^»* |x = 0)) = e-W4 /^g-iW/2t 

Y 27rt 

= -i {y\ e-'^"^ \x = 0), t > 0, (22) 

and from (!2T1) we find the result 

(y|e-*^*|x = 0) = z((?/|e-'^*|x = 0)), < t < T. (23) 

This is the main result allowing us to relate quantities calculated in the 
definite and indefinite representations. It states that, for an initial time 
interval, certain transition functions calculated in the definite and indefinite 
representations coincide up to a phase factor. The upper limit of coincidence 
T is the leftmost point on the positive time axis where the transition function 
becomes singular (a proper distribution not representable by an ordinary 
function of y). 

We now discuss the structure of the transition functions {y \ e~^^* |x = 0) 
and {{y \ e~^^^ | x = 0)) as functions of t extended to the complex plane. The 
above results do not depend on any bounds on the real part of the spectrum 
of H in either representation. As noted in the introduction to this section, 
however, the most interesting case occurs when the real part of the spectrum 
in the indefinite representation is bounded below (or if not, the spectral 
density decreases sufficiently fast) so that 

((l/|e-^^*|x = 0)) 

can be analytically extended to the lower half plane in t. We will discuss this 
case@ 

^^The discussion for other permutations of conditions on the definite or indefinite spectra 
would proceed similarly. 
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Since we are not assuming any conditions on the spectrum of H in the 
definite representation, the quantity 

(^|e-^*|x = 0), 

is not necessarily the boundary value of an analytic function in t on the range 
of t for which it exist s0 However, since 

2((l/|e-^*|x = 0)) = (y|e-*^* |a: = 0), 

on an initial interval < t < T, we may recover by analytic continuation the 
full 

((y|e-^M^ = 0)) 

for all t in the lower half of the complex plane from the values of {y \ e~*^* |x = 0) 
on this initial interval. 

In specific cases, the spectrum of H may, in addition, be bounded above in 
the definite representation. This is true for the free phantom particle and the 
phantom harmonic oscillator, but not of the two-particle system described 
in section [TOl In this case, the function 

{y\ e"*^* \x = 0) 

can be analytically extended to the upper half plane in t. Together with 

z((2/|e-^*|x = 0)), 

it therefore defines a single analytic function on the complex plane, with pos- 
sible singularities on the real axis. The value for real t of {y \ e"*-'^* \x = 0) 
is obtained by approaching the real axis from above, while the value of 
i {{y\ e~*^* |x = 0)) , is obtained by approaching the real axis from below. 
These singularities may, among other possibilities, take the form of branch 
points and cuts, as some of the examples in this article illustrate. 

There will always be at least one singularity on the real axis at the origin 
t = 0, where the transition function approaches the delta distribution. Re- 
peating the analysis of the differential equation for negative times, and using 
the fact that the free phantom particle transition function for t < requires 
the opposite sign on the right of (12^ . we see that i {{y\ e~*^* |x = 0)) will 

For an example where this function is not the boundary vahie of an analytic function, 
see section [TOl 
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instead coincide with — {y\ e~*^* \x = 0) for small negative t. As a result, in 
the special case where there are cuts on the real axis, the origin will be a 
branch point with the corresponding branch cut extending to the left. 

In general, there will be a singularity in t wherever the solution becomes a 
proper distribution not representable by an ordinary function. For example, 
the harmonic oscillator, due to the common periodicity of all the modes, will 
have singularities in t at all multiples of the half period. In this case, the 
upper limit T of the range of validity of 

{y\ e-'^' \x = {])=i{{y\ e~'^' | x = 0)) , < t < T. (24) 

is one half of the period. The fact that there is such an upper limit appears to 
be a special consequence of the even spacing of the harmonic oscillator energy 
levels, which causes the wave function to be refocused to a delta function after 
a half period. In general, perturbations that are not linear or harmonic are 
expected to destroy this even spacing, erasing the singularity that forces us to 
take T finite. Therefore, for generic potentials, except at exceptional values 
of the constants in the potential, the period may be expected to become 
infinite, so that T will be infinite and the above result should be valid for all 
positive t. 



6 Ground state expectation values 

In this section we will argue that ground state expectation values appropri- 
ate to the indefinite representation can be extracted from the corresponding 
calculation in the definite representation. 

Since a theory is completely determined by its ground state expectation 
values, one may conclude that the definite representation encodes all the 
physical information relevant to the indefinite representation, so that we can 
use whichever representation is most convenient to perform calculations!^ 

The definite representation may be preferred for functional integral calculations in 
field theories, since the functional integral obtained from the definite representation is in 
many cases manifestly covariant (modulo boundary conditions), whereas the functional 
integral obtained from the indefinite representation is not manifestly covariant due to an 
unusual integration region for the phantom components of the fields |3J [H]. Examples 
where this happens include the time-like component of a gauge field [5] and the Dirac 
boson. Dirac bosons appear in Pauli-Villars regularizations of field theories containing 
fermions and in Lee- Wick field theories, and are discussed in section [TTJ The results of 
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As in sectional we will discuss the case that is of most interest in phantom 
and ghost theories, where the spectrum of H, which is real in the definite 
representation, may be unbounded above and below. In the indefinite repre- 
sentation, the spectrum may be complex [28] . In the cases of interest to us, 
the real part of the indefinite representation spectrum will be bounded below, 
which we will assume. We also assume that the indefinite representation has 
a unique state, called the ground state, for which the real part of the energy 
spectrum attains this lowest bound. 

We will now argue that, under these assumptions, either representation 
may be used to calculate the ground state expectation values relevant to the 
indefinite representation. 

General expectation values may be generated by adding a time-dependent 
source term to the action. The argument is then based on a generalization 
of the result of section [5l whose arguments and conclusions remain valid, 
with slight modifications, if the Hamiltonian contains time-dependent linear 
source terms of the form 

IJi{t)p + \(t) g, 

where n and A are real-valued functions of t. We already know that the rele- 
vant test spaces admit time-independent linear terms of this type. However, 
the argument determining the choice of test space was based on the space 
being closed under time evolution in the limit of vanishingly small time incre- 
ments, and therefore remains equally valid when the source terms are allowed 
to depend on time. 

The further modifications needed follows: Since H in flTIJl) may 

now depend on time, we replace 

in fin?]) and subsequently. The correspondingly modified conclusion ([21]) of 
section O may now be written as 

i{{y\U{tf,U)\x = Q)) = {y\U{tf,ti)\x = ^), -T<U<tf<T. (25) 

By taking functional derivatives with respect to and A(t), we may gen- 
erate expectation values of insertions of powers of p and q in the two repre- 
sentations. These will then coincide up to a factor i if all insertions are in 
the region —T < t < T where fl2^ is valid. 

this section may be seen as a justification for the covariant treatment of these fields in the 
functional integral. 
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For simplicity, we specialize our argument to the two-point function of q. 
Similar arguments can be made for other correlation functions. The above 
result implies that 

{{y,T\Tq{tf)q{U)\x = Q,-T)) 

= (2/1 ^-mr-tf) ^^-mtf-U) q^-mu+T) |^ ^ ^26) 

provided we choose T small enough and provided that 

-T <U<tf < T. 

First, we would like to get rid of the restriction on the range of U and t/ to 
obtain a useful relation also outside the interval {—T,T). Second, we would 
like to extract the ground-state expectation value from either representation. 

Although we do not know of a general argument to get rid of the restric- 
tion on the range of ti and t/, we may do so in specific cases by one of at 
least two methods. 

First, as we remarked in section[5l for generic non-harmonic potentials one 
expects T to be infinite. In this case, the restriction on ti and tf disappears 
and we are done. 

Second, in the cases where T cannot be chosen infinite, the expressions 

((yl e-^^(^-*/) ge-*^(*/-*») ge-^^(*»+^) \x = 0)) , 

and 

-I {y\ e-^^(^-*/) ge-*^(*/-*») ge-*^(*>+^) \x = 0) , 

generally do define functions of t/ and ti also outside the range {—T,T). 
Although now these expressions have singularities on the real axis at finite 
values of T, in various cases of interest they do describe analytic functions 
of t/ and ti on the entire real line0 Indeed, there are often separate consid- 
erations from which analyticity of the two-point functions may be inferred. 
Since we have shown that these expressions coincide for sufficiently small tf 
and ti, in the case where they define analytic functions they will coincide on 
the entire region of analyticity in ti and tf. 

^^See for example the phantom harmonic osciUator of section [51 
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When a ground state exists in the indefinite representation and has over- 
lap with the state |x = 0)), the ground state expectation value may be ex- 
tracted in the standard way by writing 

{{x = o| e-*^(^-*/) ge-^^(*/-*') ge-^^(*»+^) \x = 0)) 

= ^ ((x = 0|m))e-^-^'"(^~*^') ((m|m)) {{m\q\n)) e'^^^^^f-^^'^ {{n\n)) x 

m,,n,k 

X {{n\q\ k)) e-^^'=(*'+^) {(k\k)) {{k\x = 0)) , 

and 

((x = 0| e-*^(2T) 1^ ^ 0)) = ^ ((x = 0|m)) g-'^-^^T) ^^^|^^^ ^^^i^, = g)) , 

m 

and taking the limit T — ioo of the analytic continuation in T of their 
quotient to the lower half plane 1^ This analytic continuation exists given 
the bound on the indefinite representation spectrum that we assumed in the 
introduction to this section^ 

Denoting the ground state by |0)), we obtain 

((0|g(t/)g(t,)|0))= hm 

T^—ioo 

(27) 

By (I24p and fl2B]) . we know that the definite and indefinite representations 
of the numerator and denominator on the right hand side coincide up to a 
common factor i for sufficiently small real T, and will therefore have the same 
analytic continuation to the lower half plane in T. We may therefore just as 
well write 

((ok(t/)g(t.)|o))= lim 

1 — > — %o 

(28) 



((x = 


^-iH{T-tf) q^-iH(tf-U) q^-iHiU+T) 


X = 0)) 


((x = 


g-jH(2T) 


X = 0)) 



(x = 0| 


^-iHiT-tf) q ^-iH(tf-U) q ^-iHiU+T) 


|x = 0) 


(x = 0| 


g-i_ff(2r) 


|x = 0) 



^^We have assumed here for notational convenience that the spectra are discrete and do 
not contain zero norm states, so that a basis can be chosen satisfying {{m\n)) =^ ±(5„ 
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The generahzations when this is not the case are straightforward. 

It is important to note that in general we have to take x = in the final state 
{{x — 0|, since the eigenstates \m)) do not necessarily belong to the subspace of analyti- 
cally continuable test functions. Since the indefinite representation configuration space is 
the imaginary axis, {{x\m)) is only guaranteed to be meaningful for a; = 0. This restriction 
on the final boundary condition may be relaxed in cases where all the eigenstates \m)) 
belong to the test space. 
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This formula expresses the ground state expectation value relevant to the 
indefinite representation in terms of quantities calculated in the definite rep- 
resentation, which may not have a ground state. 

To make the connection with the path integral, we note that it is only 
when — T < ti,tf < T that the expressions have a simple path integral 
representation. With this provision, we may express the result as 



We call the reader's attention to the fact that this formula enables us to 
extract the indefinite ground state expectation value from the path integral 
obtained from the definite representation, where there is no ground state. 

Once again, we repeat that for sufficiently generic potentials we expect 
T to be infinite, so that there is no need to restrict the ranges of and 
tf. When T is finite, the path integral calculation may be done assuming 
— T < ti,tf < T, and the general two-point-function for ti and tf outside this 
range must in general be obtained from separate analyticity considerations 
in ti and tf. 

Even if the upper limit for T is finite, the restriction on T may still be 
overcome even without such analyticity considerations if the definite spec- 
trum is bounded above, as in the case of the phantom harmonic oscillator 
discussed below. When there is such a bound, the expectation values on 
the right hand sides of fl27|) and fl28l) describe the boundary values of an- 
alytic functions defined respectively on the lower and upper half planes in 
T. Since these coincide for small enough T, they in fact describe a single 
analytic function with singularities or cuts on the real axis, but which can be 
analytically continued from the upper to the lower half plane past the above 
interval of coincidence. This analytic function, and its limit for T —ioo, 
can therefore be recovered from the value of the definite representation for 
any range of T. 

After this rather technical discussion, it may be helpful to remind the 
reader why the result fl29l) is useful to us. As we remarked, in concrete 
applications the definite representation is often the one where covariance is 
manifest, a property that becomes obscured in the indefinite representation. 
If we do this, the absence of a ground state in the definite representation may 



mTq{tf)q{t,)\0)) 



lim 



•too 





(29) 
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make a Euclidean time quantization problematic, whereas the real-time quan- 
tization considered in this paper remains well-defined. In a non-perturbative 
path integral approach, we do not know the ground state a priori, but the 
fixed boundary conditions in quantities such as 

{x = 0,T\Tq{tf)q{ti) \x = 0,-T) 

occurring in fl29l) are naturally and easily imposed. 

We note in conclusion that, since in both representations we have to 
perform an analytic continuation to extract the ground state expectation 
values, the indefinite representation provides no relative advantage compared 
to the definite representation for the calculation of these quantities. 

7 Inverting the construction 

So far we have argued that, for certain classes of Hamiltonians, ground state 
expectation values in an indefinite inner product representation can be ob- 
tained from transition amplitudes calculated in the positive-definite Hilbert 
space quantization of the Hamiltonian. 

It is conceivable that, for some Hamiltonians that do not belong to the 
classes we discussed, the transition amplitudes in the positive-definite quan- 
tization on a given initial time interval may be analytically extensible to the 
lower half plane. In this case, one would still expect this analytic continuation 
to display properties of the indefinite inner product representation. 

Studying this would require developing the argument of the previous sec- 
tions in reverse, and is beyond the scope of this article. 

8 The phantom harmonic oscillator 

In this section we study a harmonic oscillator with opposite sign action. It 
is now the indefinite representation that has positive energies and a ground 
state. We show that ground state expectation values for the indefinite repre- 
sentation may be extracted from the definite-representation path integral in 
the usual way, leading to an unambiguous iO prescription for the two-point 
function. 
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Let us first review the canonical quantization of the phantom Lagrangian 



X . 



(30) 



We read off that 

p — —mx, 

and canonical quantization instructs us to take [x,p\ — i. Defining as usual 



— p ] , 

2 V muj 



X p ) , 

muj 



(31) 
(32) 



where by convention 
we obtain 
and the Hamiltonian is given by 



a; > 0, 
[a, a^] = 1 



'2m 



mio 



■ X 



Note that we would have obtained the same Hamiltonian if we had started 
from the first-order form 



£ = - (a^d — a)cL) + uj a) a, 



(33) 



of the original Lagrangian. 

We now have two choices of representation. The first is the familiar 
positive inner product representation generated from a state |0) satisfying 



The resulting states 



|n) 



a |0) = 0. 

1 



(at)" |0) 



have positive inner products and negative energies 



1 



E„, = -a; I n + - ) , n = 0, 1, 2, 
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and there is no ground state. In position space, the normahzed states are 

{x\n) = (—Y^' H^{^x)e--^''\ (34) 

The other representation hves in the indefinite inner product space intro- 
duced in section [3l and is generated from a state |0)) satisfying 

at|0)) = 0, ((0|0)) = 1. 

For notational convenience, we define 

h = a\ 

so that 

[b,b^ = -l, b\0)) = 0. 

An easy consequence of the negative sign in the commutation relation is that 
the inner products of the states 

\n)) ^ ^ {b^r m 

V nl 

ahernate in sign, and that these states have positive energy 

En = +u(^n+^, n = 0,1,2,..., 

so that the state |0)) is in this ground state. From the differential 

equation implied by b\0)) = 0, we may construct these states in position 
space. We find that they are 

{{zx\n)) = -== (—) i/„(vW i-x)) e-'--- /^ (35) 

on the imaginary axis where the wave functions take their values according 
to the discussion of section El These are normalized states with respect to 
the indefinite inner product ((• • -I • • •)) defined in section [3l In fact, we have 

{{m\n)) = {-r6^n. 
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This indeed coincides with the inner product imphed by the commutation 
relations of b and 6^ The completeness relation is 



1 = Xl'^)) ^^^1^)) 

n=0 

oo 



n=0 



We can straightforwardly evaluate the ground state expectation value 

{{Tb{t)b^{0))) = -e{t)e-'''\ (36) 



Expressing 
this implies 



(a + at) = -=(6t + 6), 



((rx(t)x(0))) = -^e-l*l (37) 



Note that the sign is opposite to that of the ordinary harmonic oscillator. 

We now consider the path integral quantization of the theory. We start 
from the Hamiltonian formulation of the path integral corresponding to the 
positive-definite inner product representation, which will enable us to care- 
fully compute overall normalization factors important for the subsequent 
discussion. We need to calculate 

J [dp] [dx] expi J dt [px — H) (38) 
[dp] \dx\ expi / dt ( pi + — H — x^ \ 

= aSo / n ^ n ^ E {^^ (^^+1 ~ + ^ + 



p- , mu?_ 2 



\t^o \ 27: At J -L.-L V V 27r At J 

X exp i E ^ 



2 2 

Xj+i — Xi \ muj 2 



At 
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where we have integrated over the Pi in the last step. The prefactor comes 
from the integration over the final Pi, while the final Xj is fixed. 

Note that the last line contains the expected discretization of the La- 
grangian C in equation fl30l) . However, it is important to note that the 
measure factors 
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are different from the ones that would have been obtained from the usual 
harmonic oscillator —C, which has a factor e~^^^^ instead of e*'^/^. 

We will be interested in calculating ground state expectation values ap- 
propriate for the indefinite inner product representation, whereas the above 
path integral is the one we would obtain from the positive- definite represen- 
tation. However, in section [6] we argued that the former can be extracted 
from the latter via the usual procedure of calculating the path integral for 
real time intervals T and then taking the T — —too limit of its analytic 
continuation, as long as we take the initial and final boundary conditions 
fixed at 

Xo = xn = 0, 

and provided these states have some overlap with the ground state. This is 
exactly the same procedure one would use to extract ground state expectation 
values, if one existed, for the positive-definite inner product representation^ 
The change of variables 

= 2^(lp sm -jr- (39) 

p 

can be checked to have Jacobian 



and we obtain 



e-/^ hm xl^li fdipe 
N~.oo V 27rT J-l J '^P 

p=i 



Nm T-f r , iNm IN 



2txT V 2 



, , miV2 / p-K uo'^T^\ 2 

X exp ( —I ) — — 1 — cos — q„ 

' 2T \ N 2N^ 



^^In other words, the T —ioo Umit automatically chooses whichever representation 
has a ground state. 
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The exponent here has sign opposite to that of the usual harmonic oscillator. 
It bears repeating that this does not cause any problems in the real time for- 
malism used here0 The general definition in section [5] gives a well-defined 
meaning to the path integral as a distribution. More concretely, the above 
integrals have an unambiguous meaning as the Fourier transforms of distri- 
butions, and can be operationally evaluated in a non-perturbative setting by 
approximating the integrand by a convenient family of test functions. More 
precisely, the Fourier transform of the integrand gives a distribution that 
is locally equivalent to a unique continuous function in a neighbourhood of 
zero. The above integral is then by definition the value of this function at 



zeroo 

According to the discussion of the section O the answer, here calculated 
for real T, is expected to be analytically extensible to a function on the lower 
half plane that may have singularities and cuts on the real axis. The path 
integral expression will give the exact result as long as we keep careful track 
of overall factors of i. 

To this end, note that for each p such that 

the corresponding integral gives a factor 




p7T Uj'^T'^ 
1 — COS — — 



N 



the absolute value being superfluous here, whereas for 



pn uj'^T'^ , 



In fact, as discussed in section [Sj the integral could even have been calculated as a 
distribution in the Euclidean formalism despite the exponentially unbounded behaviour 
of the integrand in that case. 

^''The integral in each variable can easily be checked to exist as an improper Riemann 
integral, though it does not have meaning as a Lebesgue integral. However, in the present 
case where the integral is multi-dimensional, an iterated Riemann integration would not 
be invariant under general changes of variables [29j. This, and the fact that in general 
the integrand in the real-time formalism is a proper distribution, forces us to define the 
integral in terms of a transform of a distribution. 
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we obtain a factor 



P7C uj'^T'^ 
1 — cos — 



One can check that fHOl) is satisfied for no values of p as ^ oo when 

< T < tt/w, 

for the single value p = 1 when 

7l/uJ <T < 27l/uj, 

for the two values p = 1 and p = 2 when 

27i/uj <T < Sn/uj, 

and so on. We obtain 



JV^oo V 27lT 



\ 



^ / pn u^T^ 

2 1-cos- 



A^2 



2'kT 



X 



Using the identity 



cos 



El] 

n) 



\ 



N-l 



n 1 



p=i 



2m (1 - cos f ) 



N~l 

no 

p=i 



— cos — I 



N 



NJ 2^-1' 



the factor in braces is seen to be unity. Furthermore, it can be shown that, 
as A^ CX3, one gets 



N-l 



n 1 



p=i 



2N^ (1 - cos f ) 



N-l . 

p=i ^ 

sina;T 
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and we obtain the final result 



{x = 0,T\x = 0,0) = zL^-Z'^J e^-/\ l , ^, . (41) 

y 2tt\ sinu;i | 

This is the physical result, valid for purely real T. We can express it as 
the boundary value when approaching the real axis from above of a function 
analytic on the whole complex plane with cuts on the real axis as 



{x = 0, T\x = 0, 0) = e^-/\ . '7T^-nV (^2) 

Here the meaning of the square root is defined by specifying the cuts to be 
on the intervals 

(2n - 1) vr/tu < T < 212^/00, n G Z. 

The endpoints of the cuts are precisely integer multiples of the half-period. 
At these times the original position state is refocused and again becomes a 
delta distribution, so that singularities are indeed expected according to the 
general analysis of section [51 

To see that the overall powers of i are correctly reproduced, note that 
as one travels horizontally from x + iy to x + 2tt + iy for positive y, the 
value of sin {x + iy) = sin x cosh y + i cos x sinh y describes a clockwise half 
circle around the origin in C, so that ^/^\I^oi^+Jy) picks up a factor —i, and 



l/ysin (x + iy) picks up a factor i, as required. 
This has the spectral representation 



27rsincj(T + zO) V V 2 



which has the operator interpretation 

{x = 0, T\x = 0, 0) = (x = 0| e-'"^ \x = 0) 

00 

= ^ (a; = 0| e-'"^ \n) {n\x = 0) 



(43) 



n=0 

00 

^e-'^"^ |(a; = 0|n)|' 

n=0 
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in terms of the states fl34l) with positive-definite inner product (■ ■ -I • ■ ■) and 
negative energies En = — f — nu. Notice that only the states for even n 
have nonzero overlap with |a; = 0), and that the path integral gives exactly 
the correct factors | {x = 0\n) p. Also note that the powers of i required for 
correspondence with pTl) are correctly reproduced due to the common factor 



The series converges in the sense of Abel summation, which is equivalent to 
having the +iO in fH3!) . 

The positive-inner-product, negative energy spectral representation (H3l) 
is analytically extensible to the upper half plane. As expected from our 
general discussion in section El this function can be analytically extended, 
crossing the real axis on an interval to the right of the origin, to the lower half 
plane, where it can be used to obtain the indefinite inner product, positive 
energy representation. This function has cuts on the real axis and is given 
by 



2tt sin uj z 

with the constraint that the cut starting at z = is taken to the left of the 
origin. 

The values of the indefinite inner product representation inner product 
for real T are given by approaching the real axis from below, where we get 
the spectral representation 



27r sincj (T — iO) 

|^e-<t) (^-^0) + 1 e-<^+2") + ■ ■ ■ ^ , (44) 



TT 

for all real T. This has the operator interpretation 

i {{x = 0, T\x = 0,0)) = i{{x = 0\ e-'^^ \x = 0)) 

oo 

= z ^ ((x = I e~'"^\n)) {{n\n)) {{n\x = 0)) 



n=0 
oo 



5^e-^"^((n|n))|((x = 0|n))r, (45) 



n=0 



38 



in terms of the states fl35l) with indefinite inner product ((• ■ "I • • ■)) and pos- 
itive energies En = ^ + nuo. Again, the — in fl44p is equivalent to Abel 
summation. We see that the path integral has exactly supplied the correct 
prefactor i to satisfy the relation 

{x = 0, T\x = 0,0) = i {{x = 0, T\x = 0, 0)) , 

described in our general analysis of section [51 

Note that the path integral gives the correct factors 

,11,9 Prruj 
\{{x = 0\0))\' = ^ — , 

and 

{{n\n)) = +1, n even, 

since only the states for even n have nonzero overlap with ((a; = |. 

We now calculate the two-point function in the path integral approach. 
We will first do the calculation for finite time T and fixed boundary conditions 
x(0) = x(T) = 0, and then show that the correct ground state expectation 
value is unambiguously extracted from the result. 

Consider the two-point function first in momentum space 

^ p=i 

p¥=p 

X exp I —I ) — — 1 — cos — o„ > 

The integrands, including the expression on the last line, are distributions. 
Again, the integrals have a rigorous and unambiguous interpretation in terms 
of Fourier transforms of distributions. Evaluating these, we obtain the mo- 
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mentum space two-point function 



I'^c.y^^) N-m^J 



2n sin a; (T + iO) Af^oo 



27r sin a; (T + iO) \mT J (pvr/T)^ — cu^ 
Transforming to position space using (!39l) . we obtain 
(x = 0,T|Tx(t2)a;(ti) |x = 0,0) 

J2 ^ ^'"^ ^ 



27r sincu (T + iO) \fnT J ^-^^ {p-n/T) 



; mu - / 2^ \ ^ COS - cos 

27r sin a; (T + iO) VmTy ^ {p^/TY-uj^ 



-e 



j7r/4 



mc<j i / — coscu (T — 1^2 — + coscu (t2 + ^ly 



27r sin u; (T + iO) 2mci; \ sincuT 

(46) 

where we have used formula 1.445 (7) to evaluate the sums explicitly. 
The sum converges as long as — is not an integer. At these values we obtain 
poles, corresponding to the refocusing of the initial delta function by the 
harmonic oscillator potential after an integral number of half-periods, which 
makes the overlap with the final state diverge. 

It is important to note that the above summations would not converge 
for imaginary T. The path integral result must be fully evaluated for real 
T before taking the analytic continuation in T to extract the ground state 
expectation values below. 

We have obtained the two-point function for real time and fixed boundary 
conditions in the positive-definite inner product representation, which has 
no ground state. As discussed in section [6l we may therefore extract the 
ground state expectation value appropriate to the indefinite inner product 
representation from the analytic continuation to the lower half plane in T as 

//niT- 1+ \ r+ Mnw r = 0, T| Tx(t2) a;(ti) |x = 0, -T) 
( Tx t2 a; ti =^lim -. --7^, — 

T^-i<x [x = 0, T\x = 0, — T) 



2muj 
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which agrees with the operator formahsm resuh (157|) in the indefinite inner 
product representation. 
From the above and 



mu / IX \ ,+ rmo / ix 



a; + — , 



X — 



2 \ UJ J y 2 \ u 

we easily obtain the desired 

1 I" rik p-ikt 

((O|T6(t)6t(o)|O)) = -0(t)e-* = - / - ' ^ . ■ 

I J zir k — UJ + t(J 

9 The pole prescription is not derivable from 
a convergence factor 

It is important to note that the two-point function (H7|) . which we have 
unambiguously calculated starting from the real-time path integral, does not 
coincide with the naive result one would obtain if one were to replace 

Lj ^ u + ie 



in the path integral expression fl38|) . as is sometimes incorrectly argued to 
be necessary to make the integral converge. In fact, such a replacement 
is not necessary since, as we have pointed out, the integral already exists 
in the distributional sense without such an ie factor. Instead, the correct 
and unambiguous result fB7|) so obtained for the two-point function can be 
written as 

1 , I 1 f dk e-'^^ 



2muj im J 2tt — fc^ + u;^ — zO' 

which is what one would obtain by inverting the quadratic part of an action 
with the opposite replacement 

UJ UJ — iO. 



Although little practical advantage would have been gained in this example 
by doing so, we could in fact have reproduced the correct two-point function 
by adding a term —iex"^ to the original action fl38l) . Adding such a term 
would make the originally bounded integrand exponentially unbounded of 
quadratic order, so that this is not a convergence factor. As explained in 
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section despite its exponential growth, such a path integral is calculable 
as a Gel'fand-Shilov distribution, and it will indeed give the correct result as 
e ^ 0. However, the fact that this is not a "convergence factor" shows that 
the usual naive convergence argument will not work here. 

This observation is closely related to the fact that the ordinary Euclidean 
rotation of the path integral fl55]) produces an integrand of quadratic expo- 
nential growth. In other words, in contrast to the usual case, the Euclidean 
path integral is worse than the real-time path integral. Again, despite this 
bad behaviour, such a Euclidean path integral can be calculated as a distribu- 
tion using the methods of section [5l and the correct results will be obtained. 

10 Combining harmonic and phantom har- 
monic oscillators 

In this section we briefly describe the two-particle system consisting of a 
harmonic oscillator and a phantom harmonic oscillator, with Hamiltonian 

p1 moj'^ 2 vi ^^"^ 2 
2m 2 ^ 2m 2 ^ 

This system is instructive in exhibiting some features not seen in the examples 
we have considered so far. In particular, in the positive definite Hilbert space 
quantization 

• The spectrum is neither bounded above nor bounded below. 

• The transition amplitude exists for all real T except for simple pole 
singularities (not cuts), but is not the boundary value of any analytic 
function. 

• Even so, the restriction of the amplitude to a suitable initial interval 
in T can be extended to an analytic function on the lower half plane, 
and indeed to the entire complex plane except for the poles on the real 
axis. 

In the indefinite representation, on the other hand 

• The spectrum has a ground state. 
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• The transition amplitude is the boundary value of an analytic function 
extensible to the lower half plane, and indeed to the entire complex 
plane except for poles on the real axis. 

Still the arguments of section [5] will hold. The transition amplitude for the 
two representations agree for initial real T, so that ground state expectation 
values appropriate to the indefinite representation can be extracted from 
either quantization. 

To obtain the transition amplitudes, we use the following known result 
for the ordinary harmonic oscillator 



{xi = 0,T\xi = 0,0) = e^*"/^ 
which, combined with the result (H3l) 

(X2 = 0,T|X2 = 0,0) = e*"/^ 



moo 



2n sinu; (T — iO) ' 



27rsinct; (T + iO) ' 

for the phantom harmonic oscillator in the positive definite representation, 
gives 

(xi = 0,X2 = 0,T\x^ = 0,X2 = 0,0) = (-)L^-/-J 

We emphasize that, even though the system does not have a ground state, 
the real-time transition amplitude is well-defined. However, due to the factor 
(— we do not obtain the boundary value of any analytic function - 
the function cannot be analytically extended to either the upper or lower 
half plane in T, which is as expected in the absence of either bound on the 
energy. We also point out that the singularities on the real line are not cuts 
now, but poles. 

On the other hand, choosing the indefinite representation for the phantom 
oscillator, while keeping the definite representation for the ordinary oscillator, 
the combined system has a ground state. Since 



^((x2 = 0,T|x2 = 0,0)) = e^'^/y 



muj 



27rsina; (T — iO) ' 



we find 



i ((xi = 0, X2 = 0, T|xi = 0, a;2 = 0, 0)) 



In sintuT 
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where 



((0,V'|</>',^')) = (0|0)'((^l^))'- 



The transition function can be analytically extended to the lower half plane, 
as expected in the presence of a lower bound on the energy. 

Consistent with the arguments of section [5l the two results agree for an 
initial real interval, concretely < T < vr/o;, which allows us to obtain the 
result for the indefinite representation from the calculation in the definite 
representation. 

11 The Dirac boson 

In this section we discuss the quantization of a Dirac field with bosonic 
statistics. This kind of field arises, for example, as regulator fields in Pauli- 
Villars and Lee- Wick regularizations of theories containing Dirac fermions. 

This field is particularly interesting, since its Euclidean functional inte- 
grand is exponentially unbounded. Despite this, both the real-time and the 
Euclidean functional integrals may be computed using the methods of section 
[5l We will do the computation in the real-time formalism. 

Relativistic field theory requires locality (causality) and existence of a 
ground state. These conditions, imposed on a bosonic Dirac field, requires 
that we choose an indefinite inner product representation for anti-particle 
modes. However, we shall see that the indefinite representation is unsuitable 
for formulating a covariant functional integral, since there is no covariant way 
of separating these modes when we are integrating over all configurations, 
including those that are not solutions to the equations of motion. 

For this reason, it will be very useful to us that, as we have argued, we 
can extract the same results from the manifestly covariant functional integral 
based on the definite representation. 

Our starting point is the Dirac action 



where we take the field (j) to be bosonic. We first note that, because the spinor 
field is bosonic, any Euclidean version of this action [33j is unbounded below, 
so that the Euclidean path integrand is exponentially unbounded. Since the 
field appears quadratically, such an unbounded integrand could be treated 




(48) 
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using the methods of section O However, it will be easier to work in the 
real-time formalism. 

We may expand the action in modes as 

= / E / (^^^^^ + ^^^^0 ~ (^^^^^ " ^^^^^)) ' ^^^^ 

where a* and bp are bosonic. This leads to the commutation relations 

[a;, ] = {2nf 6\p - q) [b;, fof ] = {27cf 6\p - q) 6^'. (50) 

with respect to which the field satisfies locality |34j- However, if we define 
the vacuum state by 

a;|o) = 6;|o) = o, 

obtaining a positive-definite state space, it is easy to see that we may cre- 
ate states with arbitrary negative energy by applying suitable anti-particle 
creation operators bpK 

We can obtain a local (causal), positive-energy representation by keeping 
the commutation relations of the modes the same but changing the ground 
state to satisfy 

a;\o)) = b;^\o)) = o. (51) 

Since this does not affect the commutation relations of the fields and 0^, 
the desired locality is preserved. However, as in section [8], the representation 
now has an indefinite inner product. This is a general feature of local fields 
with the "wrong" statistics]^ discussed, for example, in . This choice 
of ground state is equivalent to interchanging what we mean by anti-particle 
creation/annihilation, and we could, if we wished, reflect that in the notation 
by calling bp ^ bp^ everywhere. 

Nothing prevents us from constructing a functional integral based on 
the indefinite representation fISUl) and flHTl) . However, this cannot be done 
covariantly. To understand this, note that in the indefinite representation 
(ISTj) . each mode generates a positive-definite factor in the state space, 
while each bp mode generates an indefinite factor. As a result, the two sets 
of modes would have be represented differently in the functional integral 
constructed from this representation. Specifically, as we saw in section [31 the 

^-'^But note that Pauli-Villars regulator fields are allowed to have the "wrong" statistics. 
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configuration space of a ghost degree of freedom is the imaginary axis, so 
the range of integration for the real and imaginary parts of would become 
the imaginary axis [ll[8|,[5], whereas this would not happen for a^. However, 
the mode expansion fj^9|) depends on a particular space-time decomposition 
and is not Lorentz- invariant. While the particle- antiparticle distinction can 
be made covariantly on the space of solutions, a covariant distinction is not 
possible in the space of arbitrary time-dependent configurations of the field 
used in the functional integral. The path integral resulting from the indefinite 
representation, while correct, would not be manifestly covariant. 

Since the definite representation treats all modes the same as far as their 
range of integration is concerned, the definite representation can be used 
to generate an ordinary functional integral, based on the usual sum over 
configurations using the Dirac action fHHl) . that does not suffer from this 
lack of covariance. Fortunately, we have learned how to relate indefinite and 
definite representations. 

To calculate, note that, due to the sign of the energy term, the action 
for each mode in is equivalent to that of the phantom harmonic 
oscillator (1551) . As we did in section [HI we may then extract the ground state 
expectation values for the physically relevant indefinite representation from 
the functional integral performed in the definite representation. We finc@ 

((0 iTbJit) 6,(0) 1 0)) = ^hm^ (^ = o>|4o,_T) 

= -9{t) e-'^"', 

where the calculation on the right was done in section [8] in the definite rep- 
resentation. From this, we obtain the correct Dirac two-point function, ap- 
propriate for the indefinite inner product theory, as 

{ (0 1 T 0(x) 0(y) 1 0)) = ^ hm^ (0 = o,T|0 = O,-T) 



where the quantities on the right hand side are calculated using the con- 
ventional, manifestly covariant, path integral constructed from the definite 



^Here bp corresponds to a in ((33|) . which corresponds to 6^ in the notation of ([36]) 
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representation@ 

This two-point function coincides with that of a fermionic Dirac particle. 
As required, the field satisfies locality. 

We note again that we could have reproduced this by replacing m — *■ 
m — ie in the original action fj48l) . obtaining an exponentially unbounded 
integrand that could still be calculated in the distributional sense according 
to the discussion of section O Note, however, that the correct replacement 
m — >■ m — is opposite to the naive convergence factor, and that our above 
real-time calculation showed such a replacement to be unnecessary. 

12 Complex energies and non-perturbative Lee- 
Wick type pole prescriptions 

We consider a simple example where an interaction gives rise to complex 
energies in the indefinite representation. For this example, we show that the 
path integral exists and provides a unique pole prescription for the relevant 
contour integrals. 

The example consists of a harmonic oscillator a and a phantom harmonic 
oscillator b satisfying 

[a,at] = l, [6,6t] = -l, 
with hermitian Hamiltonian 

H = LU {a^a - b^b) + 7 {a^b + b^a). 
in an indefinite inner product representation based on a state |0)), where 

a|0)) = 6|0)) = 0, ((0|0)) = 1. 

This system may be solved in the operator formalism by defining operators 

= + a = a — ib 

bi^ = a) — ib\ b = a + ib, 

that have commutation relations 

[a, a^] = 0, [6, = 
[a,b^]=2, [b,a^=2. 



^•^We should qualify this by noting that while the action used in the functional integrand 
is covariant, the boundary conditions are not. 
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Here and create null states from the vacuum |0)). The Hamiltonian 
becomes ^ ^ 

H = -{uj - i'^) a)h + -{uj + i^y) b^a, 

so that the energy eigenvalues are 

Em,n = (u; - Z7) m + (u; + ^7) n, m, n = 0, 1, 2, . . . 

The spectrum includes complex eigenvalues. Since the Hamiltonian is her- 
mitian, these come in complex conjugate pairs, as expected in an indefinite 
inner product space, and the corresponding eigenstates are null. 

We now apply the considerations of section [5] to construct the path inte- 
gral representation of this system. Since all terms in the action are at most 
quadratic in the position and momentum of either oscillator, and since no 
terms of the form Xipi or X2P2 appear, we may obtain the ground state ex- 
pectation values appropriate to the indefinite inner product representation 
from a path integral calculated as a distribution either on the test space 

corresponding to the indefinite representation used above, or on the test 
space 

corresponding to a positive-definite inner product representation for the b- 
oscillator. It was the content of section [6] that both representations encode 
the same ground state expectation values. 

As we previously discussed, the definite representation gives a simpler 
path integral. With this choice, the action used in the path integral 

= ^ J dt [a* {idt — UJ + i'j)b + b* [idt — u; — 27) a} . 

is real-valued, since the path integration measure is the usual one. 

On the other hand, in the indefinite representation, we saw that the 
configuration space for each hermitian phantom degree of freedom is the 
imaginary axis. As a result, the integration path in the path integral measure 
for each of 6 + 6* and i{b — b*) is rotated to the imaginary axis, causing the 
above action to be complex-valued. The integrand is therefore exponentially 
unbounded, in this case of linear order of growth. 
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However, despite this, it was the content of section that the the corre- 
sponding path integrals still exist as distributions on the chosen test spaces. 

Because they lacked a non-perturbative definition of the indefinite inner 
product theories corresponding to such unbounded path integrands, various 
authors [H [21 13] studied the consistency of ad hoc prescriptions for defin- 
ing the two-point functions, and more complex diagrams, order by order in 
perturbation theory. Lee and Wick [H, [2] introduced one such prescription, 
requiring the contour integrations defining various diagrams to be continu- 
ously deformed to avoid the movement of poles in the complex plane as we 
increase parameters such as 7 starting from zero. Cutkosky, Landshoff, Olive 
and Polkinghorne [5] attempted to generalize this prescription to coalescing 
singularities not covered by Lee and Wick. Boulware and Gross |1| discussed 
the problem from the path integral point of view but did not succeed in 
defining these theories non-perturbatively, mainly due to the unboundedness 
of the path integrand. 

The framework of this article allows a non-perturbative functional integral 
to be calculated in the distributional sense for various kinds of unbounded 
path integrands, including the ones studied by Lee and Wick. In principle, we 
expect the ad hoc prescriptions of the above authors to be either reproduced 
or corrected in our non-perturbative framework. 

We start with the simplest such calculation, where we derive a prescrip- 
tion of Lee- Wick type for the two-point function from the path integral. 

Consider the problem of obtaining the two-point function 



from the above action in the path-integral approach. If we were to naively 
invert the relevant quadratic term, we would obtain 



A glance at the operator representation confirms that this is wrong. The 
correct result is reproduced by 



where the integration contour T has been displaced from the real line to pass 
above the pole at 1 + i'j. This prescription, in the spirit of Lee and Wick, 



((Ta(t)6t(0))) 




ikt 




ikt 



(52) 
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seems ad hoc from a naive path integral point of view, but can be rigorously 
derived in our formalism. 

We will now show that the correct result can be unambiguously obtained, 
without any ad hoc prescription, from a careful evaluation of the path inte- 
gral. We already discussed in section [9] how the unusual pole prescription in 
the free two-point function of b follows unambiguously from the path integra- 
tion. The current discussion generalizes this to the case of complex energies. 
As we have seen, we need to evaluate the path integral for fixed boundary 
conditions between real times T and — T, and then take the limit T — icx3 
of the analytic continuation of the result to complex T. Since we have al- 
ready done this for the free oscillator in section [HI we will start from those 
results and evaluate the path integral by treating the quadratic interaction 
perturbatively in 7 starting from the form 



S 



dt (a^ {idt -uj)a-b^ {idt - cu) 6 + 7 {a^b + b^a)) 



of the action. Nothing is lost in this approach, since the perturbation is 
quadratic, so that the series can be exactly summed, giving the exact result 
for the full path integral. 

Our starting point is the result fH^ . describing the two-point function of 
the free 6-field given the fixed boundary conditions. The two-point function 
for the a-field is identical except for an opposite overall sign. Performing 
perturbation theory with fixed boundary conditions at real —T and T, an- 
alytically continuing the result in T and taking the limit as T — — 00 as 
discussed, we find 
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where the integral contour T has been deformed from the real axis into the 
upper half plane to pass above the pole at ui + i'-y to reproduce the exact 
path integral result on the previous line. There is no need to postulate this 
contour prescription ad hoc, since we have obtained it from an exact path 
integral calculation using the non-perturbatively defined path integralo 
The expectation values 

{{Tb{t)¥{0))), {{Tb{t)aHO))), {{Ta{t)¥{0))), 

may be similarly computed and put together to reproduce the exact operator 
formalism result flS21) . 

This method generalizes to more complex diagrams in interacting theories 
extending the above quadratic action. Consider, for example, an a-b loop 
diagram, expected to be proportional to 

dk f 1 1 \ / 1 1 



7 27r\/i; — Ci; + Z7 k — u ~ i'y J \p — k — uj — p — k — u + i-y ^ 

where the integration contour is so far unspecified. As above, we now do an 
exact, non-perturbative path integral calculation that determines the appro- 
priate integration contour. The path integral may be unambiguously calcu- 
lated for real T and the expectation value extracted by continuing T -—>■ —ioo 
is unique. Taking this limit term by term in the perturbation series around 
A = gives 

dk 1 1 2 C dk 1 1 



-7 



27rk — uj + iOp — k — uj + iO J 2tt {k — uj + iO)^ p — k — uj + iO 
2 r dk 1 1 
~^ J 2^k-uj + t0{p-k-uj + i0)3 ^ " ' ' 

where the integrals are over the real line. Again, we may not arbitrarily 
exchange sums and integrations. The integrals may be calculated by closing 
the contours in either the upper or lower half plane. Choosing the latter, 
only the poles contributed by the left factors contribute, and denoting 

^(^) = T^— ^' 

p — k — UJ + lO 



^^Note that we would have obtained a wrong answer if we had interchanged the order 
if integration and summation, first formaUy summing the geometric series in momentum 
space and only afterwards taking the Fourier transform. We are prevented from doing this 
by the fact that the momentum-space geometric series does not converge for the full range 
of k over which we integrate. The above answer is unambiguous. 
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the result is proportional to 

[/(a;)-7W)+yr'M----] 

- 7^ [fH - 7^ ifn^) + 1' {f)""H -■■■] 

+ 7' [/'M - f (/')"('^) + 7' - • • •] 

+ ••• 

(X f{uj + 27) + /(a; - Z7) - 7^ [/^(cu + ^7) + f{uj - i^)] +■■■ 
^ f{^ + il) /(^ - ^7) 

1 + 72 /2(c^ + i7) ^ 1 + 72 /2(^; _ 

2 1 1 

oc 1 1 

p — 2a; p — 2 (a; + i7) p — 2 (a; — ^7) ' 

since 17^/2(0; ± ^7)! < 1 so that the geometric series converge. The poles 
are precisely at the positions one would expect from the various possible 
combinations of two intermediate particles with energies uj ± i^. This result 
is reproduced by the contour integral 

r dk f I ^ \( ^ ^ \ 

Jy' 2n \k — UJ + k — uj — J \p — k — uj — p — k — uj + i^J^ 

where the integration path T stretches from —00 to 00, but is deformed to 
pass below the two poles at p — a; ± ^7 and above the two poles at a; ± i'j. 
This coincides with Lee and Wick's prescription that the contour should be 
deformed so that no poles cross it as we change 7 continuously starting from 
7 = 0. We have derived it unambiguously as a non-perturbative result in the 
path integral formulation. 

We have seen that Lee- Wick type contour prescriptions can be obtained 
unambiguously from a non-perturbative path integral. Although we have 
only considered two very simple toy examples, there is every reason to expect 
that similar non-perturbative calculations can be used to either reproduce or 
correct the prescriptions of Cutkosky, Landshoff, Olive and Polkinghorne 
unambiguously. 

13 Conclusion 

In this article, we investigated the non-perturbative quantization of phantom 
and ghost degrees of freedom by relating their representations in definite and 
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indefinite inner product spaces. For a large class of potentials, we argued 
that the same physical information can be extracted from either representa- 
tion, and provided a non-perturbative definition of the path integral for these 
theories. We applied the results to the study of ghost fields of Pauli-Villars 
and Lee- Wick type, calculating non-perturbatively previously ad hoc pre- 
scriptions for Feynman diagram contour integrals in the presence of complex 
energies. 

The initial motivation for this work was to understand the nonpertur- 
bative path integral quantization of certain bosonic Pauli-Villars ghosts for 
which the Euclidean path integrand is unbounded [36] • These ghosts were 
an important ingredient in the construction of a class of functional integral 
measures that are generally covariant, background- independent and confor- 
mally invariant. The current article provides the missing ingredient to make 
the construction of ^36j rigorous. 

The mathematical framework of this article may have applications in the 
study of phantoms appearing in certain cosmological models of dark energy. 
It may also be of interest in the approaches to the cosmological constant 
problem based on phantoms, or in the approaches based on symmetries con- 
sisting of rotations of the configuration space to the imaginary axis. 

We pointed out some pitfalls in trying to derive ie prescriptions in Feyn- 
man integrals from naive convergence terms in the action. We showed that 
for ghost fields, the correct prescription is in fact opposite to what one would 
obtain from such convergence factors. 

Finally, the methods introduced here may conceivably be useful in the 
complex-coordinate or complex-momentum approaches to the study of reso- 
nant scattering. 

Acknowledgments 

We would like to thank Prof. Antal Jevicki and the Brown University Physics 
department for their support. 

References 

[1] T.D. Lee and G.C. Wick, Negative metric and the unitarity of the S- 
matrix, Nuclear Physics B9 (1969) 209-243. 



53 



[2] T.D. Lee, A finite theory of quantum electrodynamics, in Elementary 
processes at high energy, part A, Ettore Majorana 1970 International 
School of Subnuclear Physics, Erice, July 1-19, Editor: A. Zichichi, Aca- 
demic Press, New York and London, 1971. 

R.E. Cutkosky, P.V. Landshoff, D.L Olive and J.C. Polkinghorne, Nu- 
clear Physics B12 (1969) 281-300. 

D.G. Boulware and D.G. Gross, Lee-Wick indefinite metric quantization: 
A functional integral approach. Nuclear Physics B233 (1984) 1-23. 

Seiji Sakoda, Euclidean path integral of the gauge field - holomorphic 
representation, arXiv.org preprint |hep-th/0501205 2005. 

LM Gel'fand and G.E. Shilov, Generalized functions: Spaces of fun- 
damental and generalized functions. Vol. 2, Academic Press, Orlando, 
1958. 

LM Gel'fand and G.E. Shilov, Generalized functions: Theory of differ- 
ential equations. Vol. 3, Academic Press, Orlando, 1967. 

H. Arisue, T. Fujiwara, T. Inoue and K. Ogawa, Generalized Schrddinger 
representation and its application to gauge field theories. Journal of 
Mathematical Physics 22 (1981) 2055-2059. 

R. R. Caldwell, A Phantom Menace? Gosmological consequences of a 
dark energy component with super-negative equation of state. Physics 
Letters B545 (2002) 23-29. 

James M. Cline and Sangyong Jeon and Guy D. Moore, The phantom 
menaced: Gonstraints on low- energy effective ghosts. Physical Review 
D70 (2004) 043543. 

Antonio De Felice, Mark Hindmarsh and Mark Trodden, Ghosts, insta- 
bilities, and superluminal propagation in modified gravity models. Jour- 
nal of Gosmology and Astroparticle Physics 08 (2006) 005. 

Jens Kujat, Robert J. Scherrer and A. A. Sen, Phantom Dark 
Energy Models with Negative Kinetic Term, arXiv.org preprint 



astro-ph/0606735, 2006. 



54 



Mariusz P. Dabrowski, Glaus Kiefer and Barbara Sandhoefer, Quantum 
phantom cosmology, Physical Review D74 (2006) 044022. 

Recai Erdem, A symmetry for vanishing cosmological constant in an 
extra dimensional toy model, Physics Letters B621 (2005) 11-17. 

Recai Erdem, A symmetry for vanishing cosmological constant: Another 
realization. Physics Letters B639 (2006) 348-353. 

Gerard 't Hooft and Stefan Nobbenhuis, Invariance under complex 
transformations, and its relevance to the cosmological constant problem. 
Classical and Quantum Gravity 23 (2006) 3819. 

David E. Kaplan and Raman Sundrum, A Symmetry for the Cosmolog- 
ical Constant, Journal of High Energy Physics 0607 (2006) 042. 

I. Antoniadis and E. Dudas and D. M. Ghilencea, Living with ghosts and 
their radiative corrections, arXiv.org preprint liep-th/0608094, 2006. 

W. Pauli and F. Villars, On the invariant regularization in relativistic 
quantum field theory. Reviews of Modern Physics 21 (1949) 434-444. 

Tore Berggren, On the use of resonant states in eigenfunction expansions 
of scattering and reaction amplitudes. Nuclear Physics A109 (1968) 265- 
287. 

Roger, G. Newton, Analytic properties of radial wave functions. Journal 
of Mathematical Physics 1 (1960) 319-347 . 

A. Buchleitner, B. Gremaud and D. Delande, Wavefunctions of atomic 
resonances. Journal of Physics B27 (1994) 2663-2679. 

Nimrod Moiseyev, Quantum theory of resonances: calculating energies, 
widths and cross-sections by complex scaling. Physics Reports 302 (1998) 
211-293. 

G. Hagen, J.S. Vaagen and M. Hjort- Jensen, The contour deformation 
method in momentum space, applied to subatomic physics. Journal of 
Physics A37 (2004) 8991-9021. 

Y.K. Ho, The method of complex coordinate rotation and its application 
to atomic collision processes. Physics Reports 99 (1983) 1-68. 



55 



[26] R. de la Madrid, Description of resonances within the rigged Hilhert 
space, arXiv.org preprint quant-ph/0607168 , 2006. 



[27] Carl M. Bender, Sebastian F. Brandt, Jun-Hua Chen, Qinghai Wang, 
Physical Review D71 (2005) 025014. 

[28] J. Bognar, Indefinite inner product spaces, Springer- Verlag Berlin, 1974; 
A.I. Mal'cev, Foundations of linear algebra, W.H. Freeman and Com- 
pany, 1963. 

[29] P. Cartier and C. DeWitt-Morette, A rigorous mathematical foundation 
of functional integration, in Functional integration: Basics and applica- 
tions, eds. C. DeWitt-Morette, P. Cartier and A. Folacci, Proceedings 
of NATO Advanced Study Institute held September 1-14, 1996, Plenum 
Press, New York, 1997. 

[30] G.W. Johnson and M.L. Lapidus, The Feynman integral and Feynman's 
operational calculus, Oxford University Press, Oxford, 2000. 

[31] G. Vitah, Rend. R. Istor. Lombardo (2), 36 (1903) 772-774. 

[32] G. Vitali, Ann. Mat. Pura Appl. (3), 10 (1904) 73. 

[33] P. van Nieuwenhuizen and A. Waldron, On Euclidean spinors and Wick 
rotations. Physics Letters B389 (1996) 29-36. 

[34] M.E. Peskin and D.V. Schroeder, An introduction to quantum field the- 
ory, Perseus Books, Cambridge MA, 1995. 

[35] N.N. Bogolubov, A. A. Logunov, A.I. Oksak and I.T. Todorov, General 
principles of quantum field theory, English edition, Kluwer Academic 
Publishers, Dordrecht, 1990. 

[36] Andre van Tonder, Worldsheet covariant path integral quantization 
of strings, to appear in International Journal of Modern Physics A, 
arXiv.org preprint hep-th/0606017^ 1996. 



56 



